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Effect of periodic potential on the stripe phase realized in the higher Landau levels is investigated 
by the Hartree-Fock approximation. The period of the potential is chosen to be two to six times of 
the fundamental period of the stripe phase. It is found that the stripe aligns perpendicularly to the 
external potential in contrast to a naive expectation and hydrodynamic theory. Charge modulation 
towards the Wigner crystallization along the stripe is essential for the present unexpected new result. 
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When the two-dimensional electrons are in a strong 
magnetic field such that the lowest Landau level is par- 
tially filled, the fractional quantum Hall effect or the com- 
posite fermion Fermi liquid state is realized at low tem- 
peratures depending on the filling factor. On the other 
hand, it is almost established experimentally |l|-|3| and 
theoretically Q-Q] that the stripe phase or the unidirec- 
tional CDW state is formed in ultra-high mobility sam- 
ples at low temperature when the third lowest or higher 
Landau levels is nearly half filled. The orientation of 
the stripe is fixed relative to the crystal axis in the per- 
pendicular field, although in-plane magnetic field can af- 
fect it. [p[-p0[ Hence, the longitudinal resistivity is much 
larger in the [1,1,0] direction than in the [1,1,0] direction. 
According to a theoretical calculation of the conductiv- 
ity tensor of the stripe phase, |11| this indicates that the 
stripe runs along the [1,1,0] direction. 

The reason why the stripe has preferred orientation 
has not been known. It is hard to believe that the crys- 
tal structure affects the orientation of the stripe. In the 
experiments the magnetic field is a few teslas, so the wave 
function of each electron in the third Landau level is 
spread over more than several hundred angstroms. So 
details of the crystal lattice structure will be averaged 
out. On the other hand, investigation of the surface mor- 
phology of high mobility samples shows structure running 
along the [1, 1, 0] direction. |l^ , |l3|] It is possible that this 
structure orients the stripe. 

To sec if this morphology is relevant or not, we consider 
effect of a periodic potential of wave vector Q p on the 
stripes of fundamental wave vector Q s in this letter. We 
calculate the ground state energy in two configurations: 
in one of the configuration Q p and Q s are parallel to each 
other and in the other they are perpendicular. Naively, 
it seems that the parallel configuration has lower energy. 
As a matter of fact, if the sizes of the wave vectors are 
the same, the lowest energy is achieved when they are 
parallel to each other. Moreover, Luttinger liquid theory, 
or hydrodynamic theory, |ll],[l4|,[l5| which is considered to 
be valid in the long wavelength, low energy limit, favors 
parallel orientation. This is because in the perpendicular 



orientation, only one type of deformation of the stripe 
is caused. It is modulation of the width of the stripes. 
On the other hand, in the parallel orientation, two types 
of the deformation of the stripe are caused at the same 
time. The additional deformation is the displacement of 
the stripes. This additional coupling to the displacement 
brings more energy gain in the parallel orientation. We 
will come back to this point later. 

Now if Q p and Q s are parallel, we are in trouble. In 
the experiments, the wave vector of potential modulation 
Q p should be in the [1,1,0] direction, and the resistivity 
is lower in the same direction, which indicates Q s is in 
the [1,1,0] direction. In order to resolve this paradox we 
present in this letter results of microscopic Hartree-Fock 
calculation of the stripe phase in the external potential 
modulation. The calculation is done for the third lowest 
Landau level at half-filling. Our results show that con- 
trary to the long wavelength theory the wave vectors of 
the stripe and the potential are orthogonal to each other 
at least when the wavelength of the potential modulation 
is a few times of the period of the stripe. 

Our Hamiltonian consists of two parts: The Coulomb 
interaction part for the third Landau level is given as 
follows: 



H ° = ^E u (9)[^(^ 2 ^ 2 )] 2 ^x-g), 



v(q) 



2eq 



X ± =X±-l 2 q y , 0<X 



2irl 2 



-i < L. 



(1) 



(2) 



(3) 



(4) 



In these equations p(q) is the density operator of the guid- 



ing center, L n (x) is the Laguerre polynomial, I — ^Jfi/eB 
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is the magnetic length, and L — > oo is the linear dimen- 
sion of the system. The other part is the interaction term 
with the external potential, 

H! = ~V eKt p(q)L 2 (~q 2 l 2 )eM^o), (5) 

q=±Q P 

where V ey ± is the strength of the potential, and r*o is the 
origin of the potential. 

In the Hartree-Fock theory |ll|[l7j we introduce order 
parameters 

A(Q) = J2( a x + a xJcxp(-iQ x X), (6) 
x 

where X± = X ± (l/2)l 2 Q y . The wave vectors of the 
order parameters are expressed as 

Q = UQ° x ,kQ° v ), (7) 

with j and k being integers. We recognize a CDW state 
as the stripe phase when the order parameters with Q = 
±Q S are dominant. 

In the present calculation we generalized the procedure 
detailed in ref. (l7) to the case of rectangular symmetry, 
Q° x ^Q y . Then when 

NQ° x Q° y l 2 = 2Mir, (8) 

with N and M being integers, the Landau level splits 
into N Hofstadter bands, and we fill half of the bands 
to get the filling factor v = 1/2. First we searched the 
optimum stripe state by changing the integers N, M and 
the ratio Q x /Qy- As a result we found that the choice 
g° = Q° = Q s = N = 6, and M = 1 gives 

almost the lowest energy stripe state, the energy being 
0.30969e 2 /47reZ per electron. In this calculation we have 
retained the order parameters with wave vector Q such 
that \Q\l < Qmax' = V200. The Yoshioka-Lee's order 
parameter sum rule, |l7j 

^ = ]T|A(Q)| 2 , (9) 

Q 

is satisfied by 99.994%, which indicates that the con- 
tribution from the discarded order parameters is negli- 
gible. The obtained stripe state is degenerated in the 
direction of the fundamental wave vector as expected. 
Namely, we obtained ground state where the fundamen- 
tal wave vector of the stripe Q s is either (Q s , 0) or (0, Q s ). 
The density along the stripe is modulated slightly with 
wave number 3Q S corresponding to the linear alignment 
of electrons. Namely, the unidirectional CDW state has 
"2/cp instability" towards the Wigner crystallization, and 
energy gap is created at the Fermi level as has been first 
pointed out by Yoshioka and Fukuyama. Jl6| However, 
the modulation is not large enough to destroy the topol- 
ogy of the stripe. p8[ 



Next we take into account the external potential. 
We consider external potential with wave vector Q p = 
(Q p ,0). We apply this potential to the stripe state that 

has Q s = (Qs, 0) or (0, Q s ), respectively. The wave num- 
ber of the potential Q p is chosen to be Q s /m, where m 
is an integer. Then we calculate the order parameters 
A(jQp,fcQ s ) with j and k integers self-consistently, and 
obtain the ground state energy in the periodic potential. 
Since Q x is decreased to Q p , N becomes 6m. Namely, 
the energy spectrum consists of 6m bands in this case. 
As m becomes larger so does both the size of the matrix 
we must diagonalize and the number of order parameters; 
it takes longer time to obtain the self-consistent solution. 
Therefore, we calculated only up to m = 6. In the course 
of the calculation all the order parameters whose wave 
vector is smaller than 10/ Z is retained except for m = 6, 
for which Q max = \/40/Z. We have checked at smaller m 
that Qmax = \/40/Z is enough: although the ground state 
energy changes as we decrease Q max , the energy gain in 
the external potential changes only slightly. At m = 2 
the difference between the Q max = V200/Z and y/iO/l 
cases is largest and about 2.5%, however, at m = 3 it is 
less than 1%. For m = 4 and 5 the difference between 
the Qmax = 10/Z and -y/40/Z cases is much smaller. 

The result for the case of Q p = Q s is trivial. Once the 
external potential becomes finite, we obtain only a self- 
consistent solution where the fundamental wave vector 
of the stripe coincides with that of the external poten- 
tial. In this case the ground state energy lowers almost 
proportional to V eX f We will not discuss this case any 
more. 

For 2 < m < 6 we have obtained two different self- 
consistent solutions. Namely, in one solution Q s is per- 
pendicular to Q p , and in the other they are parallel. The 
energy gain depends on the origin of the external poten- 
tial, so optimized position is chosen. In all the cases the 
energy gain is almost proportional to V 2 xt , as long as V ex t 
is smaller than the Hartree-Fock self-consistent potential 
of the CDW, ~ 0.3e 2 /47reZ. We express the difference as 
follows: 

8E ee E(V cxt ) - E(0) = -\eV^ + 0(V c l t ), (10) 

where E(V cx t) is the ground state energy per electron 
in V^xt- The proportionality constant e is given in Ta- 
ble. 1. As shown there the most stable configuration is 
that where the wave vectors are perpendicular to each 
other. This result is opposite to that of the long wave- 
length theory. 
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FIG. 1. The real space CDW pattern in the external po- 
tential. In this figure Q p = Q B /3 and T4 x t = 0.2e 2 /Aivel. The 
Fourier transform of the (p(q)) is shown for an area of linear 
dimension of four times of the stripe period A = 2n/Q s . The 
brighter region shows higher density, (a) for the case of paral- 
lel configuration, and (b) for the perpendicular configuration. 

The pattern of the stripe in the potential is shown in 
Fig.l for the case of m = 3. When Q p _L Q s , there 
is only width modulation. |L9| On the other hand, the 
displacement of the stripe is the main deformation in 
the parallel case. For different choice of the origin of the 
potential, the width of each stripe also changes. Actually, 
there is the electron-hole symmetry, so the displacement 
of the electron stripe means width change of the hole- 
stripe, and these two patterns give the same energy. 

In order to extrapolate our result to longer wavelength, 
we calculated what we call effective bulk modulus k de- 
fined by the following equation for both configurations: 



£(Qp) = 



2irl 2 



v{Q p ) + 2irl 2 k' 



(11) 



This k is plotted in Fig. 2. This plot clearly shows that 
even in the limit of Q p — > 0, k for the parallel configu- 
ration is larger than that of the perpendicular configura- 
tion, and the energy gain is larger in the latter case. 
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FIG. 2. The effective bulk modulus k is plotted versus Q p l. 
The dashed lines are guide to the eyes. The open and closed 
circles are for the parallel and perpendicular configurations, 
respectively. 



Now we explain how the hydrodynamic theory gives 
the opposite result. When the stripe runs parallel to 
the i-axis, or Q s is in the y-direction, the hydrodynamic 
theory describes the deformation of the stripe by the fol- 
lowing Hamiltonian. n,0,na 



H = 



2L 2 



Kqi)\ux\ 



+ [«(?) + X -1 ]l<KjH 2 + 2im e C ? j / u^n_ 9 -}, (12) 

where, u* gives the displacement of the stripe, Sn^ is the 
coarse-grained component of the electron density modu- 
lation due to the width change of the stripe, \ is the com- 
pressibility from which the effect of the Coulomb inter- 
action is removed, v(q) is the Coulomb potential, eq.((|), 
m is the electron mass, and Y, K, C are parameters. 
The compressibility is related to the ground state energy 
per electron, E, by the following equation: 



2ttZ : 



dv 2 ' 



(13) 



If q is small we can equate 5n^ and p(q)L2(q 2 l 2 /2). By 
adding eq.(|5|) to this Hamiltonian, we obtain the energy 
gain in the presence of the external potential as follows: 
When the wave vector of the potential is Q p = (Q p , 0) _L 

Q s , 
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Because of the coupling between the displacement and 
width modulation of the stripes, the latter has lower en- 
ergy. J2(| This is opposite to our result. The compress- 
ibility is estimated from the ground state energy to be 
X~ 1 /2nl 2 ~ — 0.32e 2 /47reL This value is quite close to 
the effective bulk modulus k for the perpendicular con- 
figuration shown in Fig. 2. Thus our result agrees with 
that of the hydrodynamic theory only in the perpendic- 
ular configuration. 

The origin of the present behavior is the charge mod- 
ulation along the stripes, which is correlated between 
stripes. To understand it, we first consider what hap- 
pens when the filling factor is changed. There are two 
orthogonal ways to change the filling factor. In one way 
the width of the stripe is changed at fixed periodicity and 
in the other the periodicity is changed at fixed width of 
the stripes. |2l]] If the filling factor is changed uniformly, 
width change at almost fixed periodicity occurs. This is 
because the stripe periodicity is determined by the mini- 
mum of the Hartree-Fock potential. Q Next we consider 
unidirectional compression of the stripes. If there is no 
charge transfer between the stripes, the compressibility 
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will depend on the direction. If they are compressed 
along the stripe direction, the periodicity is kept con- 
stant, and the width of each stripe changes. This is easy 
direction, and the compressibility is high. On the other 
hand, in the perpendicular direction the period changes 
and the width is kept constant. This is hard direction, 
and compressibility is low. In the presence of charge 
transfer between stripes, this kind of compression will 
never be realized. 

However, the situation is different if periodic potential 
modulation of finite wavelength is applied to the present 
stripe phase, which is actually a rectangular CDW state 
with charge modulation along the stripes. The charge 
modulation shows the tendency toward the Wigner crys- 
tallization, and the period is inversely proportional to the 
stripe width. Inter-stripe correlation favors the stripes to 
have the same width in the direction perpendicular to the 
stripe. In the perpendicular configuration, Q p _L Q s , this 
correlation does not interfere the stripe to respond to 
V^ x t: the stripe changes the width as shown in Fig. 1(b). 
Thus the bulk modulus in eq. ( |TT| ) is approximately given 
by x _1 /27r/ 2 . On the other hand, in the parallel config- 
uration, Q p || Q s , stripes try to respond to V ey ± keeping 
the width the same as neighbors. So the main response is 
the displacement of their positions, and the periodicity is 
modulated. The energy gain is smaller in this situation. 
Actually, the bulk modulus of the period change is esti- 
mated to be k ~ le 2 /47rd, which is consistent with Fig. 2. 
The reason that k decreases as Q p decreases is that slight 
width modulation is allowed at longer wavelength. 

The reason why m = 2 case has low k even for the 
parallel configuration is explained by this scenario. In 
this case the displacement of the electron-stripe means 
the width change of the hole-stripe at fixed periodicity. 
This is unique to m — 2, since for m > 3 the displacement 
of the electron-stripes causes both the displacement and 
width change of the hole-stripes. Therefore, only m = 2 
has low bulk modulus in the parallel configuration. 

From the argument above we propose to modify the 
hydro dynamic theory such that the compressibility x has 
q dependence. Another fact ignored in the hydro dynamic 
theory is that the displacement of the stripe induces local 
charge modulation as pointed out above. Namely, density 
modulation Srig — iq y riQU^ is induced by u^, where no is 
mean electron density. Hydrodynamic theory including 
these effects should be constructed. 

Ishikawa and Maeda also calculated the effect of the 
external potential. |||| They applied the second order 
perturbation theory to their Hartree-Fock ground state. 
Their conclusion is that there is no energy gain in the par- 
allel configuration, and the perpendicular configuration 
has lower energy. Although the qualitative conclusion is 
the same as ours, we doubt validity of their calculation. 
As we have shown here the coupling to the external po- 
tential always deforms the stripe state, and lowers the 
energy of the system. The lack of the energy gain in 
their parallel case indicates that the intermediate states 



are not correctly taken into account. Thus their conclu- 
sion that the perpendicular configuration is favored has 
no foundation. 

In conclusion we have shown that the stripes orient 
themselves perpendicular to the long wavelength external 
modulation to gain energy without loosing correlation 
intrinsic to the stripe phase. 

Part of the numerical calculation is performed at the 
ISSP, the Univ. of Tokyo. This work is supported by 
Grant-in- Aid No. 12640308 from MEXT, Japan. 
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TABLE I. The value of ee 2 /47reZ, which is the coefficient 
of the energy gain per electron in the external potential. The 
size of the wave vector of the external potential is Q p — Q s /m, 
the direction of which is either parallel or perpendicular to the 
stripe's wave vector. 



m Q p l parallel perpendicular 

2 0.5117 0.568 0.577 

3 0.3411 0.252 0.372 

4 0.2558 0.208 0.273 

5 0.2047 0.175 0.217 

6 0.1706 0.152 0.180 



5 



